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SIGNIFICANCE  AND  EXPLANATION 


OM  THE  DOMAZH  SPACE  FOR  COMSTITOTIVE 
LAWS  IN  LINEAR  VISCOELASTICITY 

Michael  Renardy 

1*  Introduction 

The  most  general  constitutive  law  for  an  Isotropic,  inccs^ressible 
linearly  viscoelastic  fluid  has  the  form 

t 

(1.1)  T(t)  «  /  a(t-s)Jf(t,s)ds, 

•mm 

where  T  denotes  the  extra  stress  tensor,  and  J  denotes  the  linearized 
relative  strain  tensor. 

A  recent  paper  by  Saut  and  Joseph  [4]  addresses  the  question  of 
specifying  suitable  daa«in  spaces  for  the  history  of  j.  The  class  of 
"admissible"  kernels  a  is  then  given  by  the  dual  of  the  domain  space. 
Depending  on  the  choice  of  domain  space,  the  kernels  may  be  functions  or,  more 
generally,  distributions.  Saut  and  Joseph  [4]  discuss  various  choices  for  the 
domain  and  characterize  the  resulting  dual  spaces. 

In  this  paper,  a  different  point  of  view  is  adopted.  Rather  than 
specifying  a  domain  space  at  the  outset,  I  shall  investigate  the  ia^lication 
of  certain  phenomenological  restrictions,  which  are  commonly  accepted  in 
rheology.  It  will  turn  out  that  these  restrictions  alone  guarantee  a  certain 
regularity  of  the  kernel  a,  thus  implying  that  the  constitutive  law  is 
defined  on  a  certain  minimal  dcxaaln  space. 
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We  adopt  the  following  two  conditioner  which  are  easily  motivated  by 
physical  intuition  (to  knowledge,  they  were  first  formulated  by  Boltzmann 
[1]): 

1.  The  fluid  resists  deformation;  positive  strains  lead  to  positive 
stresses. 

2.  The  strain  resulting  from  a  more  recent  instant  of  time  always  has  a 
greater  Influence  than  that  resulting  from  a  more  distant  time. 

If  a  were  a  smooth  function,  these  conditions  «K>uld  be  written  as 
a>0,a'<0.  We  adopt  a  generalized  form  of  this,  which  is  made  precise 
below. 

We  shall  find  that  the  only  admissible  non- function  contribution  to  the 
kernel  is  that  which  corresponds  to  the  Me%;tonian  fluid.  In  particular, 
fluids  of  hi^er  grade  are  excluded  (they  do  not  satisfy  condition  1). 

In  chapter  3  I  investigate  the  stability  of  fluids  of  higher  grade. 
Expanding  an  argument  of  Joseph  [3] ,  I  find  that  no  fluid  of  grade  higher  than 
four  can  have  a  stable  rest  state.  Moreover,  if  the  domain  occupied  by  the 
fluid  is  infinite,  the  initial  value  problem  for  such  fluids  is  ill-posed. 
Therefore,  even  for  smooth  initial  histories,  such  fluids  cannot  be  used  to 
approximate  unsteady  flows. 
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2.  Characterigatlon  of  Admissible  Kamels 


A  priori  It  will  only  ba  aasumad  that  the  constitutive  law  Is  defined  on 
the  following  space  of  test  functions: 

B  >  e  C  [0,  «•),  >app  is  cosset,  ^(0)  ■  0}. 

m 

Thus  the  mapping  'P  J  a(s>^(s)ds  Is  a  continuous  linear  functional  on  B. 

0 

The  conditions  1  and  2  above  are  reformulated  In  the  following  distributional 
form: 

m 

1*.  /  a(s)^(s)da  >  0  for  ^  e  B  such  that  ^  >  0. 

0 


2*.  /  a(s)^'(s)d8  >  0  for  ^  e  B  such  that  ^  >  0. 

0  supp  p  C  (0,»)  compact. 

These  conditions  alone  guarantee  that  the  kernel  a  can  actually  be  defined 
on  a  much  larger  domain  space  than  B. 

Theorem: 

If  1*,  2*  hold,  then  a  can  be  extended  to  a  linear  functional  defined 
on  the  space 

F  -  e  L^(0,<*),  ^(0)  ■  0  and  <P  Is  differentiable  at  O). 

This  extension  Is  continuous  In  the  sense  that  for  any  a  >  0  there  Is  a 
constant  c(a)  such  that 

m  m  La  I 

l/  a(sV  (8)ds|  <  C(o)  (  /  1^  (s)  Ida  sup  +  1^'  (0)  |  } 

0  0  se(0,al  ” 
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(If  ^  e  F,  then  the  right  hand  side  of  this  inequality  is  finite  for  a 
email  enough)* 

Proof; 

Let  US  first  consider  how  the  kernel  a  acts  on  test  functions  f 
supported  on  an  interval  [a, ci  >  0.  For  any  such  Interval  there  is  a 
function  e  E  such  that  >  0  and  -  1  on  the  interval  [  a#  fi] . 

We  thus  have 


•  max  |^(t)|  <  w(8)  < 
teta,ei 


*  max  |w(t)|. 
tela,  M 


From  condition  1*,  it  follows  that 


1/ a(s)^(s)ds|  <  max  |w(t)|  •  /  a(s) *  g(s)ds, 

0  teto,fiJ  0 


whenever  s\q>p  <fi  is  contained  in  la, 0}.  Hence,  the  kernel  a  can  be 
extended  to  a  continuous  linear  functional  on  the  space  I  a,  0] .  According 
to  a  well  known  result  [2],  this  is^lies  that  a(8)ds  ~  duls),  where  M  is  a 
finite  regular  countably  additive  measure  on  [a, 0].  Because  of  condition 
2*,  a*  (the  generalised  derivative  of  a)  has  the  same  property.  This 
laq>lies  that  a  is  a  function  of  bounded  variation.  It  is  easy  to  conclude 
from  1*  that  a  is  positive.  Moreover,  we  have  -a'(s)d8  ■  dlj(s)  for  a 
certain  regular  countably  additive  measure  w.  The  characteristic  function 
Xj  of  any  open  interval  z  C  (a, 0)  can  be  approximated  by  a  sequence  of 
positive  test  functions.  Property  2*  and  Lebesgue's  theorem  guarantee  that 
M(I)  is  positive.  Hence  u  is  a  positive  measure,  which  implies  that  a 
is  monotone  decreasing.  (We  have  a(t)  -  a(t')  ■  -M(t,t')  for  almost  every 
t,t'). 


I'  T.  -f. — ^  ■ 


I 


I 

i 

I 


I 


Let  us  now  consider  test  functions 

be  such  that  4*  (s)  “  s  for  s  6  [0,a] 
a 


'P  supported  on  [0,a]. 

and  4*  >  0«  Then 

o 


Let  4'5,  e  E 


-4>  max  K'(t)|  <  <  4*  max  |^'(t)|. 

®  tetO,a)  ®  te[0,a] 

a  a 

Therefore,  we  have  /  a(s)<^(s)ds  /  b(s)^'(s)ds,  and  b(s)ds  >  dv(8)  for  a 

0  0 

finite,  regular,  countably  additive  measure  v.  since  w  already  know  how  the 
kernel  a  acts  on  functions  supported  away  from  0,  we  need  only  be 
concerned  with  the  possibility  that  v  has  an  atom  at  0.  We  thus  obtain 

«D  «» 

(2.1)  /  a(s)^(s)ds  >  C  ^'(0)  +  /  i(s)^(s)ds  for  ^  6  B. 

0  0 

where  S  is  a  positive,  monotone  decreasing  function  defined  on  (0,*). 

Condition  implies  that  C  >  0,  and  for  the  integral  to  be  defined  on  B 

it  is  necessary  and  sufficient  that  i(s)*s  is  locally  integrable.  (Since 
0» 

”  e  C  ,  ^(0)  s  0,  we  have  in  general  a  as  s  >  0).  From  here  the 

theorem  is  immediate. 

Remark; 

When  speaking  of  "domain  spaces",  we  have  regarded  (1.1)  as  a  functional, 
i.e.  as  a  mapping  of  histories  of  J  to  present  values  of  T.  k  different 
sometimes  useful  view  of  convolution  integrals  such  as  (1.1)  is  that  of  an 
operator,  i.e.  a  mapping  of  histories  of  J  to  histories  of  T.  With  this 
Interpretation,  the  domain  space  for  the  relation  (2.1)  will  be  different 
from  F  (of  course,  it  depends  on  the  space  that  the  history  of  T  is 
required  to  be  in). 
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Remark  on  the  formulation  of  constitutive  eauationa; 

1  have  shown  that,  1£  Boltzmann's  theory  is  reformulated  in  a 
distributional  setting,  then  the  Newtonian  fluid  appears  as  a  special  case. 
Constitutive  models  that  combine  memory  Integrals  and  Ne%^onian  contributions 
(or  nonlinear  modifications  thereof)  are  therefore  nattural.  Such  models  are 
usually  used  for  polymer  solutions  with  Newtonian  solvents  and  have  also  been 
suggested  for  polymer  melts.  On  the  other  hand,  fluids  of  higher  grade,  when 
taken  as  "exact”  models  rather  than  singular  perturbation  expansions,  are 
Incompatible  with  Boltzmann's  postulates. 
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3.  Monexlatence  of  "materials  of  the  differential  type*  [S] 


As  we  have  seen,  fluids  of  grade  higher  than  1  are  not  consistent  with 
Boltsmann's  conditions.  In  this  chapter  1  show  that,  in  addition,  fluids  of 
grade  higher  than  4  can  not  be  stable  at  rest  if  the  domain  is  sufficiently 
large.  X  conclude  that  such  materials  can  not  exist.  For  orders  two,  three, 
and  four,  stability  of  the  rest  state  is  possible  only  if  the  first  normal 
stress  difference  has  a  sign  incompatible  with  experimental  values.  These 
results  extend  those  of  Joseph  [31  who  showed  that  a  fluid  of  grade  n  can  be 
stable  only  if  the  coefficients  of  ^  and  have  the  same  sign. 

As  shown  in  [31 ,  the  stability  problem  for  an  nth  grade  fluid  at  rest 
leads  to  the  characteristic  equation. 


(3.1) 


o  •  A  k(  o) 

A 

-  vK  *  «,»♦  —  ♦ 


where  A  are  the  eigenvalues  of  the  stokes  operator  for  the  domain  occupied 
by  the  fluid  and  k(0)  is  a  polynomial  of  degree  n  -  1.  The  A  can  take 

A 

any  positive  value  if  the  domain  is  chosen  appropriately. 

Consider  the  limit  A  ■*'0  (which  occurs  if  the  domain  grows  large). 

A 

Then  one  root  of  (3.1)  tends  to  sero  and  is,  in  a  first  approximations,  given 

by  0  -  A  .  The  other  n  -*  2  roots  tend  to  •*  and  their  first 
0  A 

approximation  is  a  ■  A_a  or  «  -r — —  .  If  n  is  greater  than 

m  n“i  no, 

m  n-1 

four,  at  least  one  of  the  roots  to  this  last  equation  will  have  positive  real 

part.  (This  argument  remains  valid  even  if  an  integral  term  as  considered  by 

Saot  and  Joseph  [4)  is  added  to  the  constitutive  law.  Such  a  term  is  only  a 

'*owe*^  ,.der  perturbation).  If  n  is  three  or  four,  we  conclude  from  the 

preceding  argument  that  stability  is  only  possible  if  at  and  o  ,  are 

U  11*1 
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Asu^aitj  jua 


negative.  According  to  [3] ,  must  then  also  be  negative.  This  proves 

my  claim  for  n  •  3.  For  n  ■  4,  let  us  consider  the  limit  A  ♦  In  this 

IR 

case  the  roots  are  approximated  by  those  of  the  equation 

(3.2)  )c(o)  -  o  +  0^0^  +  ■  0. 

We  already  know  that  a^,  are  negative.  There  Is  at  least  one 

negative  root  whence 

k(o)  -  (o  +  o  )(e  +  +  B,o^) 

(3.3)  012  3 

-  BlOg  +  <®2  • 

We  can  conclude  that  8^  and  must  be  negative.  It  Is  then,  however, 

easy  to  see  that  the  quadratic  polynomial  has  a  root  of  positive  real  part 

unless  B^  also  negative.  As  a  consequence,  »  B^  +  B^Oq  Is  negative. 

It  must  be  noted  that,  for  n  >  4,  we  not  only  find  eigenvalues  of 

positive  real  part,  but  for  A  >0  the  real  part  tends  to  Infinity.  In 

in 

Infinite  domains,  the  Initial  value  problem  for  such  fluids  Is  therefore  Ill- 
posed,  l.e.  solutions  generally  do  not  exist  even  for  smooth  Initial  data. 
Rlvlln-Erlcksen  fluids  of  grade  higher  than  4  therefore  are  not  only 
Inadmissible  as  "exact"  models,  but  also  as  approximate  models  for  unsteady 
flows,  they  are  an  unstable  approximation. 

The  phenomenon  that  an  approximation  to  the  equation  does  not  always 
yield  an  approximation  to  the  solution  Is  well  known  In  numerical  analysis, 
however.  It  Is  by  no  means  limited  to  numerical  methods.  Any  kind  of 
approximation  must  be  assessed  for  Its  stability.  It  Is  my  conclusion  that 
there  are  two  poss'  jle  lines  for  further  research t 
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1 .  Abandon  the  Idea  of  retarded  motion  expansions  [6]  and  seek 


different  and  hopefully  stable  approximation  schemes. 

2.  Rlvlln-Erlcksen  type  approximations  might  be  "saved"  If  additional 
approximations  are  Introduced.  This  situation  Is  quite  common  In 
numerical  analysis.  Many  discretization  methods  require  that 
spatial  and  temporal  mesh  sizes  satisfy  certain  Inequalities.  This 
means  In  particular  that  the  two  approximations  must  not  be  made 
Independently . 


I 
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